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Abstract. By means of a modified hypervirial theorem we derive simple expressions 
for the integrals of products of Airy functions. Present results contain earlier ones as 
particular cases. 

Airy functions appear in many problems of mathematical physics. For that reason 
some time ago Gordon [1] and Albright [2] derived integrals of products of Airy functions 
that do not appear in commonly available tables of integrals and functions. However, 
they are necessary for several problems of physical interest [1,2]. 

The purpose of this letter is to propose a simple approach for the derivation of such 
integrals. It is based on a straightforward modification of the well known hypervirial 
theorems [3] and may also be useful for other similar problems. 

Let A(x) and B(x) be solutions to the eigenvalue equations 

LA = aA, LB = bB 

L =D 2 -x, D = 4- (1) 
dx 

We can easily prove that for any linear operator O 

J A[L, 0}Bdx=(a-b) J AOB dx + W(A, OB) (2) 

where [L, O] = LO — OL is the commutator, W(u, v) = uv' — u'v is the Wronskian and 
the prime indicates differentiation with respect to x. 
We need the following two cases 

[LJD] = f"D + 2f\L + x) + f (3) 

and 

[L,g]=g" + 2g'D (4) 

where / and g are differentiable functions of x. It follows from equation ([2]) with O = ft) 
and O = g that 

J Af'DB dx + 2bj Af'B dx + 2j Axf'B dx + J AfB dx = (5) 
(a - b) J Af'DB dx + W(A, fB') (6) 
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and 

J Ag'DB dx = J AgB dx - ^ J Ag"B dx + ^W(A, gB) (7) 

respectively. Notice that we can use the latter equation twice with g = f and with 
g' = f to remove the two integrals in the former equation that contain DB = B' . Thus 
we derive the following expression 

J AfB dx + 2 J Axf'B dx - - J Af'B dx+(a + b) J Af'B dx 

- / AFB dx = W(A, fB') - lw(A, f'B) 

+ ^W(A,FB) (8) 

where F = f f dx. From this master equation we easily obtain any integral of the form 
/ x n A(x)B(x) dx such as the particular cases considered by Gordon [1]. Notice that he 
chose functions of the form A[a(R + j3i)] but we can set the scale a = 1 without loss of 
generality so that (3\ = a and 02 = b in our notation. For example, if F = 1 we obtain 
his equation (B.10): 

f A(x)B(x) dx = -—[A'(x)B(x) - A{x)B'(x)} (9) 

J (X 

If we choose a = b = then the master equation (jBJ) gives us all the integrals of the 

form I n {yi,y2) = f x n yiy2dx considered by Albright [2]. For example, when / = x n we 

obtain the recurrence relation 

n,/ Nr,/ \ , n(n-l)(n-2) r n 8 , „ . . , x n+1 A(x)S(x) 
x n A(x)S(x) da; = — — - — - — - / x n ~ 3 A(x)B(x) dx + 



2(2n + l) J v / v / 2n + l 

n(n - l)x n ~ 2 A(x)B(x) nx 11 ' 1 [A(x)B'(x) + A'(x)B(x)} 
2(2n+ 1) + 2(2n+ 1) 

_ x"A!(x)B(x) n = t (10) 
2n + 1 

On the other hand, equation ([7]) enables us to connect the integrals I n (yi,y' 2 ) with the 
-fm (2/1,2/2) ones. Finally, if we take into account that 

we easily express the integrals hiy'iiU^) m terms of those discussed above. 
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In our opinion the modified hypervirial theorems provide a most convenient simple 

and systematic approach to the calculation of integrals of products of Airy functions. 
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